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: ABSTRACT 

' The shape of a galaxy is constrained both by mechanisms of formation (dis- 

sipational vs. dissipationless) and by the available orbit families (the shape 
^ ' and amount of regular and stochastic orbits). It is shown that, despite the 

tH- ' often very flattened shapes of banana and fish orbits, these boxlet orbits gen- 
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0^ ' constraint from the shape of fish orbits is relaxed at (large) radii where the 
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erally do not fit a triaxial galaxy in detail because, similar to loop orbits, they 
spend too little time at the major axis of the model density distribution. This 
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c/3 . 1 INTRODUCTION 



density profile of a galaxy is steep. 



Destruction of box orbits and growth of chaos in triaxial potentials with a divergent central 
force are the classical arguments for adopting axial symmetry in modeling galactic nuclei 
(Binney & Gerhard 1985). The role of the so-called boxlet orbits is comparatively less un- 
derstood theoretically. These centro-phobic orbits, named after various food families such 
as bananas, fishes and pretzels (Miralda-Escude & Schwarzschild 1989), are flattened in a 
way grossly similar to that of the box orbits in a cored potential, but for reasons not en- 
tirely clear they appear not as useful as box orbits when it comes to fit a flattened galaxy 
model. Pfenniger & de Zeeuw (1989) suspected this is because a closed banana orbit comes 
in only a limited range of axis ratios. Schwarzschild (1993) suggests the reason being that 
the boxlets have their density maxima at the turning points (corners) rather than on the 
major axis. Similar qualitative arguments include that a thick banana orbit and its mirror 
image together have twice as many corners as a box orbit, which makes them less flexible in 
fitting a smooth density model (Kuijken 1993). Syer & Zhao (1998) studied banana orbits 
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in the separable non-axisyinmetric disc models of Sridliar & Touma (1997), and found that 
the self-consistency is first broken down near the symmetry axes of the model. 

The present analysis is also motivated by the question: what are the main factors in 
preventing a galaxy being triaxial? Commonly one might attribute to (a) a dissipational 
formation of a rotation-supported galaxy, or (b) the growth of stochasticity and shrinking 
of available phase space for low resonance boxlet orbits in triaxial models once a steep cusp 
or a central blackhole is developed. If these were the only factors, one would expect that 
if any strongly triaxial models exist, they should be the ones with zero stochasticity and a 
shallow cusp. Yet Syer & Zhao found that the continuous family of banana orbits in the fully 
separable discs of Sridhar & Touma cannot make a self-consistent flattened model (density 
flattening of about 0.6) even when the surface density profile is as shallow as /?~°°^. On the 
other hand Schwarzschild (1993) and Merritt (1997) studied models with a E6 flattening and 
a much steeper R~^ projected cusp and found that they are nearly self-consistent despite 
large stochasticity in the phase space in non-separable models. There is clearly a third factor 
in determining the shape of a model, which is (c) the shape of low resonance boxlet orbits. 
The mismatch of the shape of banana orbits excludes them as suitable orbits for building a 
triaxial model even when there are unrestricted number of them (Syer & Zhao 1998). Same 
can be said about loop orbits. A fish orbit, with one of its apocenters aligned with the major 
axis, is arguably ideal for building a self-consistent model, but somewhat puzzlingly they 
are populated merely at the level of 20% in Schwarzschild (1993) models. Their role has not 
been studied in recent literature as detailedly as the banana orbits and the stochastic orbits. 
Currently about a few dozen grid-based numerical models have been built to explore the 
parameter space of triaxial models, while a few thousand might be adequate, however, this 
might be just on the edge of present numerical power. 

On the analytical front, a recent study by Zhao, Carollo & de Zeeuw (1999) shows that the 
shape of orbits very near the major axis of a potential is very useful for classifying whether 
an orbit is helpful or not helpful for building a potential. They solved for the equations 
of motion infinitely close to the major axis analytically. They found that boxlet orbits are 
often ill-suited for building the potential because of mismatch of their shape or curvature 
at the major axis with that of the density model; they are not fiat enough. This is rather 
surprising because these boxlets are often more flattened than that of the model, both in 
the globally averaged sense. Another surprising result is that the problem of mismatch of 
curvature of boxlet orbits is more severe for shallow- cusp ed (but not cored) systems than for 
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steep-cusped systems; this goes differently from tlie expected trend from numerical studies 
if the stochasticity is the main factor in limiting galaxy shapes. However, their analysis is 
rigorous only for 2-dimensional scale-free non-axisymmetric discs. 

Here we seek to highlight the generic problem of boxlet orbits at the major axis (§1.1). 
We then limit ourselves to fish orbits, but extend the analytical method of Zhao et al. to 
a general triaxial potential (§2). This allows us to explore the full range of triaxial models 
of galactic nuclei and judge whether fish orbits are well-suited or ill-suited for building a 
specific model. Finally we make the connection of our results with properties of observed 
galaxies and with previous works (§3) and present our conclusions (§4). 

1.1 General problem of boxlets and loops: pendulum dilemma 

The major axis of a model is always where the amplitude of the angular momentum of a 
regular orbit J{t) reach a local maximum. This is determined by the equations of motion: the 
torque from the flattened potential always increases the angular momentum when the orbit 
comes close to the major axis, and decreases it when leaving the major axis. The situation is 
analogous to a pendulum oscillating around the major axis (see Fig. 0), which never spends 
enough time near the desired position (the major axis), and always too much time near the 
undesired positions (the turn-around points). Hence while stars need to slow down so as to 
spend more time near the major axis to fit the model, the equations of motion tell them to 
do the opposite. This is perhaps the main limiting factor when boxlet orbits and loop orbits 
are used in making models (Zhao et al. 1999). 

The best that we can do to prolong the stay of a star near the major axis is to let the star 
reach its apocenter (a local maximum in radius) at the major axis, for example, as in the 
case of a fish orbit (Fig. |l]) or the like such as the 3 : 4 pretzels. The situation is analogous 
to a pendulum with a variable length, e.g., a spring attached by a weight becomes longer 
at the equilibrium point, and moves slower in terms of the angular speed. Nevertheless the 
following arguments show that the effect is still limited. To be specific we concentrate on the 
fish family, but the same arguments apply equal well to any boxlet family with an apocenter 
on the major axis. 
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2 CAN FISHES GET OUT OF THE PENDULUM DILEMMA? 

Consider a general triaxial model with a density p{r,9,ilj) and potential $(r, ^,-0) in a 
spherical coordinate (r, 9, i/j) where the short axis (z-axis) is the pole with ^ = 0, and the 
azimuthal angle ifj runs from in the yz plane to | in the xz-p\a.ne. Define the parameters 
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so that they describe the fiattening of the potential in the xz and xy planes at radius tq on 
the major axis (x-axis with 6' = -0 = |); generally e^ > e^ > since the contours are fiatter 
in the xz plane than the xy plane. 

Now launch a fish orbit tangentially from its major-axis apocenter at radius rg with an 
angle (pitch angle) x from the xy plane and an instantaneous angular momentum Jq- Let 
the time t = at the launching, then the trajectory can be approximated as 
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for a short enough period of time with t -^T = —, where T is the local dynamical time scale 
and f^ = rdr^{r, 9, 'i/;)\r=ro is the local circular velocity squared. The convex dependence of 
time for both J{t) and r(t), meaning that t = is a maximum, is a result of the torque from 
the fiattened potential and our choice of the apocenter. 

As in the Schwarzschild (1979) method we compute the amount of time that an orbit 
spends in a set of cells and compare it with the amount of mass in the cell as prescribed by the 
volume density model. First we divide the density model into shells with equal logarithmic 
intervals radially [rj = Alogr). Then tessellate the angular dimension into equal solid angles 
with a size AQ = /3^, where f3 is the typical angular scale of the cell. Hence the amount of 
mass in the cell, 

dlogp 



AM = p{r, 9, 'i/jyr]AQ oc r'^~"'''S{9, i/j), 1 + 7^ 



dlogr 



(5) 



where the volume density p(r, 9, i/j) is approximated as a power-law of slope 1+7^ near radius 
r, and the angular shape of the volume density is specified by a shape function S{9, ip). We 
shall pretend that 7^ is a constant in the subsequent derivation (as in a scale- free model), 
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but generalize it to be a shallow function of r at the end; in the same spirit we treat the 
flattening parameters eg and e^ of the potential model. 

Now our orbit crosses the cells with an angular velocity uj{t) so the amount of time At 
that it typically spends in the cell with a characteristic angular size jS is given by 

A^=-, ^(^) = ^' (6) 

where ^ is a geometrical constant of order unity and J{t) is the instantaneous value of the 
amplitude of the total angular momentum of the orbit when crossing the cell. Here we have 
ignored the possibility that our orbit may cross a cell radially by arguing that a fish orbit 
takes much less time to cross a cell tangentially than radially. 

Comparing At with the mass in the cell AM, we find that the ratio of the two for an 
orbit [r{t),9{t),ip{t)\ is given by 
At a{t) _r''^it) 



AM 5(^(t),VX^))' J{t) ' 

where a{t) can be taken as the effective projected density of a fish orbit in the direction 

[9{t),ilj{t)]. The important thing here is that while S{6{t), '0(t)) is generally at maximum on 

the major axis for any realistic density model, the projected density of a fish orbit a{t) can 

actually be at minimum at the same place if 



— = [ee cos X + e.^ sm xj - 7a. I 1 - IsT:^ ) (8) 



'qVc, 



is positive or zero, where we have substituted in eqs. @ and (^. In other words a fish orbit 
will be a "bad" orbit if its angular momentum is high with 
4 > ^ _ Max[e^,6e] 

Eq. (^ is in fact generally applicable for spotting "bad" orbits. For example, it predicts 
that loop orbits are universally "bad" since they always visit the major axis with an angular 
momentum Jq slightly greater than that of a circular orbit, tqVc- 

Our result on fish orbits can also be casted to an even simpler form: eq. (^) is always 
satisfied with a positive l.h.s. and a negative r.h.s. if e^ > 7^, i.e., the iso-potential contours 
is flatter than what fish orbits can support. This region is marked as "bad- fish-zone" in 
Fig. 0. 
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3 IMPLICATIONS TO OBSERVED GALAXIES AND COMPARISON 
WITH PREVIOUS WORKS 

Where are observed galaxies in Fig. 0? One tricky point of connecting our result with obser- 
vation is that the orientation and intrinsic semi-axes of their three principal axes (assuming 
a triaxial model) are often unobservable. Nevertheless we can set useful limits since pro- 
jected contours can only be rounder than intrinsic and their projected central cusp slope is 
generally equal to or slightly steeper than the cusp-slope 7^ as rigorously defined in eq. @, 
i.e., deprojection would shift the fish-like symbols in the scatter diagram down and to the 
right (roughly the direction pointed by the fish symbols), deeper into the "bad-fish zone". 
About half of the galactic nuclei sample are in this zone where fish orbits are "bad" for 
making a self-consistent triaxial model. 

Another tricky point is that the density profiles of real galaxies steepen towards large 
radii. Giant ellipticals, for example, have a volume density power-law slope 1 -|- 7^ in the 
range of — 1.3 in the nucleus, and 2 — 4 at radii comparable to one effective radius. 
Early studies (e.g., Binney 1978) of the kinematics and isophotes of these systems argue 
convincingly that they are generally triaxial objects, at least outside the nucleus, supported 
by radial anisotropy from box or boxlet orbits (Schwarzschild 1979). If we let 7^ vary with 
radius, our result (eq. ^ would suggest that a fish orbit is more effective in providing the 
anisotropic support to the potential at large radii when the density profile is steeper and 
the iso-potential contour is rounder than at small radii. 

Ryden (1999) in a recent preprint on statistical study of isophotes of a similar galaxy 
sample came to the conclusion that while isophotes at the inner core of giant ellipticals are 
consistent with them being randomly oriented oblate objects, the shapes become non-oblate 
at sufficiently large radii where the mixing of stochastic orbits is incomplete over a Hubble 
time (Merritt & Valluri 1996). It is reassuring to see both stochastic orbits and fish orbits 
follow the same general trend, i.e., they seem to be more helpful for supporting a triaxial 
potential at large radii than at small radii. 

Previous numerical studies of scale-free models by Schwarzschild (1993) and double- 
power-law models by Merritt (1997) with an inner volume density cusp 1 + 7^ = 2 show 
that boxlet orbits cannot support an ellipsoidal model fiatter than E6 or E7. We find that 
these models are indeed too fiat for fish orbits (cf. Fig. |^). 

It is less clear whether boxlets become more helpful or less helpful for triaxial models with 
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a shallow cusp. While Merritt & Fridman (1997) suggest that shallow cusp triaxial systems 
are easier to build than steep cusp systems by extrapolating the findings of two ellipsoidal 
double-power-law models (with 1 + 7^ = 1 or 2 at the center), our result (eq. §) predicts that 
fish orbits become "worse" when the density profile becomes shallow (1 + 7^^ < 1). Similarly, 
while Schwarzschild (1979) finds that triaxial models with a core and a finite force at the 
center can be built with centro-philic box orbits, Syer & Zhao (1998) find that separable 
non-axisymmetric disc models with an E5 shape in potential cannot be made with banana 
orbits even if the projected density profile is as shallow as R~^'^^. Using a similar analysis, 
Zhao et al. (1999) also found that the curvatures of boxlet orbits at the major axis could not 
reproduce that of a 2-dimensional scale-free non-axisymmetric disc if the power-law slope of 
the disc is small but non-zero. They suggest the existence of a forbidden-zone, very similar to 
our "bad-fish-zone". Whether the above findings are contradicting remains to be studied. It 
depends on whether fish orbits or banana orbits are essential for building a triaxial galaxy. 
If they occupy only a very small fraction of the phase space, then one might argue that 
they are marginally interesting orbits and has little effect on the shape of the galaxy. It also 
depends on whether it is crucial to include additional constraints coming from fitting density 
off the major axis and from the amount of phase space allocated to regular orbits vs. that 
allocated to stochastic orbits; these constraints are not exploited here merely because they 
are not as easy to tap analytically as the conditions of fish orbits near the major axis. 

4 CONCLUSIONS 

In summary we explore the role of fish orbits analytically in a parameter space much larger 
than possible with previous techniques. Our technique is extended on the basis of the scale- 
free disc models of Zhao et al. (1999) to a general potential. We find that fish orbits are 
generally not suitable for building triaxial models because their tendency to "speed" at the 
major axis often dominates any gain from putting them at a apocenter on the major axis. 
In essence centro-phobic fish orbits are just as "bad" as loop orbits in the sense that neither 
can spend much time at the major axis, which excludes them as helpful orbits for a triaxial 
model. Since all boxlets and loops suffer the so-called pendulum dilemma (cf. § |1.1| ), with fish 
orbits being among the most promising of all boxlet orbits to soften the pendulum dilemma, 
absence of "good" fishes for some very fiattened models imply that we have no "good" orbits 
left for making these models self-consistent. 
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The general impression from our result on fish orbits and previous works is that a shallow 
cusp and absence of chaos do not necessarily allow a triaxial model, and triaxiality is difficult 
to support for steep cusp systems as well as for shallow cusp systems. If all galaxies have 
central black holes, then most of the available phase space in a strongly flattened, triaxial 
potential will be taken by stochastic orbits, loop orbits, and high resonance boxlet orbits. Ba- 
nana orbits, fish orbits and other low resonance food orbits can come only in limited fiavors, 
i.e., axis ratios and sizes. The lack of a continuous varieties of shapes in the building blocks 
could be a serious limitation for building a smooth self-consistent models since any model 
built would likely to show undesirable spiky features reminiscent of the sharp boundaries 
of individual banana or fish. The mismatch of shapes must be another important limiting 
factor for triaxial galaxies, in addition to dissipational formation and growth of chaos. 

HSZ thanks Tim de Zeeuw for comments on an early version of the paper which helped 
greatly to improve the presentation here. 
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Figure 1. The upper panel shows two types of centro-phobic orbits: a loop orbit (thick filled band) and a fish orbit (dots) 
together with the iso-density and iso-potential contours (heavy solid and dashed lines respectively). The lower panel shows the 
radius r(t) and angular momentum \J{t)\ of a star vs. the shape function S{r{t), 0{t), ip{t)) = p{r(t), 9{t), 'il){t))/ p(r{t), 7r/2, tt/2) 
along the same two orbits, where each dot is one time step of the orbit and Sit) — > 1 as the orbit approaches the major axis. A 
pendulum with a variable length is also sketched in the upper panel. Note that \J(t)\ peaks, hence the star spends least time, 
at the major axis for boxlets, like it does for loops and the pendulum. They are thus ill suited for building a fiattened model 
unless the model density along the orbit p((r(t), S(i), -0(4)) is in fact at a minimum on the major axis which could happen if 
p(r, 0,-0) falls steeply with radius and the major axis is an apocenter. 
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Figure 2. The parameter space for scale-free ellipsoids: the cusp-slope 7^ versus the flattening 1 — c/a. Overplotted (fish-like 

symbols) are the projected central cusp slope vs. ellipticity of 56 Space Telescope observed galaxies taken simply from the 

"typical" values given in tabulated form in Faber et al. (1997). Three dashed curves, in descending order, are for triaxiality 

2_, 2 n 

parameter °^_ ^ = [0.9,0.5,0.1]; they mark the boundary of the "bad-fish zone" (cf. eq. H). Galaxies marked by the fatter fish 

symbols (those near or to the right of the dashed lines) are unlikely to host fish orbits which spend enough time on the major 
axis; their names are listed in the order of cusp slope and ellipticity with ones to the left right being "cored" and rounder than 
E2. Galaxies to the left of the vertical (dotted) line are very round models; the chance for a highly triaxial and highly flattened 
galaxy to project into the range left of the vertical line is fairly small; the probability is less than 3% for a galaxy with intrinsic 
axis ratio (1 : 0.75 : 0.5) to appear rounder than El and only 10% to appear rounder than E2 in projection from a random 
orientation. Between the dotted line and the dashed lines is a grey zone where it is more difficult to predict the effect of the 
fishes. 
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